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ANALYSIS OF TURBULENT FLOW AND HEAT TRANSFER IN NONCIRCULAR PASSAGES 
By Robert G. Deissler and Maynard F. Taylor 


SUMMARY 

Previous work on turbulent beat transfer and flow in tubes was gen- 
eralized and applied to flow in noncircular passages of equilateral tri- 
angular and square cross section. Expressions for eddy diffusivity that 
bad been verified for flow and beat transfer in tubes were assumed to 
apply in general along lines normal to a wall. Velocity distributions, 
wall sbear-stress distributions, and friction factors, as well as wall 
beat -transfer distributions, wall temperature distributions, and average 
beat-transfer coefficients, were calculated. In addition, results from 
a previous analysis for axial flow between rods were compared witb new 
experimental data. For calculating wall temperature distributions, uni- 
form beat generation in tbe passage wall and uniform beat transfer at 
tbe outer surface were assumed. Tbe application of tbe results is re- 
stricted to moderately small peripheral rail temperature variations. 
Calculations were made for Reynolds numbers from 20,000 to 900,000 and 
Prandtl numbers from 0.73 to 300. 

Results show that velocities, shear stresses, and beat transfer in 
the region near the corner were lower than average values and were zero 
at tbe corner. Friction factors and average Nusselt numbers were lower 
than in a tube. 


INTRODUCTION 

Until recently, most analyses for flow and beat transfer in passages 
were confined to circular tubes or to parallel plates. Those passages 
were analyzed extensively because of their importance in technical ap- 
plications and because their simplicity makes them amenable to analysis. 

In recent years tbe problems associated witb the use of noncircular 
passages in beat exchangers have become important. In reference 1, wall 
temperature distributions for turbulent flow in rectangular and triangular 
ducts were calculated by using experimental velocity distributions and 
average beat-transfer coefficients, together witb an assumed similarity 
of tbe wall beat-transfer and wall sbear-stress variations; no attempt was 
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ma de to calculate either the heat-transfer coefficients or the velocity 
and temperature distributions in the fluid field. Some calculations of 
velocity and shear-stress distributions in corners are reported in 
reference 2. 

As a part of an investigation being conducted at the NACA Lewis 
laboratory on heat transfer and flow in passages of various shapes, 
turbulent axial flow between rods and in eccentric annuli was analyzed 
in references 3 and 4. In the present investigation, the calculations 
were carried out for flow in passages having square and equilateral tri- 
angular cross sections. As in the previous analyses, expressions for 
eddy diffusivity that had been verified for flow and heat transfer in 
tubes were assumed to apply in general along lines normal to a wall. 

This assumption, although appearing to be reasonable , should of course 
be checked experimentally. Also, the peripheral temperature variation 
was assumed to be small compared with the variation of temperature across 
the passage. 


BASIC EQUATIONS AND ASSUMPTIONS 

Thft differential equations for fully developed flow shear stress 
and heat transfer can be written in the following form: 

t = (H + pe) ~ (l) 

q = - (k + pc p e h ) ^ (2) 

where e and % are the eddy diffusivities for momentum and heat 
transfer, respectively, the values for which depend on the amount and 
kind of turbulent mixing at a point. (All symbols are defined in appendix 
A.) In these equations, y is taken as the perpendicular distance from 
the wall, and t and q are measured in planes parallel to the wall. 
Equations (l) and (2) can be considered as definitions of e and e^. 

They can be written in dimensionless form as: 

+ — 1 -\— ( 5 ) 

Po P-o/PQ / dy + 



*0 \ k o ^0 Po c p 0 


€ \ dt + 
*575/ dy + 


( 4 ) 


where the subscript 0 refers to values at a wall and the properties in 
the definitions of u+, y+, and t+ are evaluated at the wall. 
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Egressions for Eddy Diffusivity 


In order to use equations (3) and (4) , the eddy diffusivity e must 
he evaluated for each part of the flow. As in reference 5, in the region 
away from the wall the mechanism for turbulent transfer is assumed depend- 
ent only on the velocities in the vicinity of the point, measured relative 
to the point, or on the space derivatives of the velocity. In the region 
close to the wall, the turbulence is assumed dependent on quantities meas- 
ured relative to the wall, that is, on u and y and on the kinematic 
viscosity p/p. On the basis of dimensional analysis and simplicity, it 
is assumed in reference 5 that, for the region close to the wall (y"*" < 26), 


— n^uy l 1 - e v J 


( 5 ) 


where the constant n has the experimentally determined value 0.124.^- 


In the region away from the wall (y + > 26), e is assumed to be de- 
pendent on the relative velocities in the neighborhood of the point. 

Prom a Taylor’s series expansion for u as a function of y and z. 


„ _ ^ ( &u d^u du S^u d^u \ 

= X s ? ^ * * •' S’ 5? • • 7 

where y and z are measured in normal directions in the cross section 
of the passage. Since, in the case of flow through a tube or between 
parallel plates, the velocity-gradient line is a line that at each point 
is normal to a constant -velocity line, the derivatives in the z direc- 
tion are zero. For flow in noncircular passages, the velocity-gradient 
lines near a wall are also normal to the surface, but they are usually 
curved in the center part of the passage. Inasmuch as the greatest 
changes of velocity with respect to distance take place in layers near 
the wall, the effect of the derivatives with respect to z will be 
neglected. It seems reasonable to expect that near the center of the 
flow passage the effect of the derivatives with respect to z would be 
to increase the turbulence and flatten the profile in that region. 
However, the normal turbulent profile (derivatives with respect to z 


^The quantity in parentheses in equation (5) becomes important only 
for heat transfer at Prandtl numbers appreciably greater than 1. For 
Prandtl numbers on the order of 1 or less, or for the calculation of 
velocity profiles, e = n^uy (the value of n differs from that in eq. (5)) 
and is a good approximation for the region close to the wall (ref. 6) . 
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absent) is already very flat in that region, so that the increased turbu- 
lence should not produce significant changes in the values of the veloc- 
ities. Therefore, the expression for e for y + > 26, obtained by using 
dimensional analysis, is 


e = x 2 (^/dy ) 4 5 

(d^u/dy^) 


( 6 ) 


where x has the experimental value 0.36 (ref. 6) ; this is von Karman’s 
expression. 


Further Assumptions 

In order to integrate equations (3) and (4), the following assump- 
tions are made, in addition to those concerning the eddy diffusivity 
(eqs. (5) and (6)): 

(1) The fluid properties are considered constant. The analysis 
could be carried out for variable properties, but the complexity would 
be increased. 

(2) The eddy diffusivities for momentum e and heat transfer 

are equal (the ratio a = l) . Previous analyses for flow in tubes based 
on this assumption yielded heat-transfer coefficients that agreed with 
experiment (refs. 5 and 7). At low Reynolds or Peclet numbers (Pe = RePr) , 
a appears to be a function of Peclet number (ref. 8); but for Reynolds 
numbers above 15,000, as in the present analysis, a is nearly constant 
for gases and for fluids with Prandtl numbers greater than 1. 

(3) Along lines normal, to a wall, the variations of the shear 
stress t and heat-transfer per unit area q have a negligible effect 
on the velocity and temperature distributions. It is shown in figure 12 
of reference 7 that the assumption of a linear variation of shear stress 
and heat transfer across a tube (t or q ■ 0 at passage center) gives 
very nearly the same velocity and temperature profiles as those obtained 
by assuming uniform shear stress and heat transfer across the passage. 

(4) The molecular shear-stress and heat-transfer terms in equations 

(3) and (4) can be neglected in the region away from the wall (y + > 26) 
(ref. 7, fig. 14). 



MCA TN 4384 


5 


Generalized Velocity and Temperature Distribution 

With the foregoing assumptions, equations (3) and (4) were integrated 
in reference 5. The results are reproduced in figures 1 and 2, These 
curves give the relations of u + , y + , and t + that are used in the fol- 
lowing calculations for flow in noncircular passages. 


CALCULATION OF VELOCITY DISTRIBUTIONS FOR FLOW 

■> 

i IN NONCIRCULAR PASSAGES 

For applying the relation between u + and y + in fi glare 1 to the 
calculation of velocity distributions for noncircular passages, an iterative 
procedure must be used, inasmuch as the lines of velocity gradient (lines 
normal to constant-velocity lines) are unknown at the outset. A typical, 
flow passage of the type considered is shown in the following sketch: 
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The first step in obtaining the velocity distribution is to draw the 
line of symmetry; the velocities on either side of this line are lower 
than the velocities on the line. In the present case, the line can be 
drawn immediately from symmetry considerations, as shown in the sketch. 

The next step is to draw assumed velocity-gradient lines. By using these 
assumed velocity -gradient lines, lines of constant velocity can be calcu- 
lated, as will be shown. A new and more accurate set of velocity-gradient 
lines can then be drawn, inasmuch as they must be normal to the constant- 
velocity lines. With a little practice, the velocity -gradient lines can 
be estimated quite accurately the first time, so that usually the iterative 
procedure need not be used more than once or twice. 


Calculation of Lines of Constant Velocity 

For calculating velocities at various points in the passage, it is 
convenient to define the following velocity parameter: 



This parameter is used in place of u + because the shear stress in the 
definition of u + varies with position. Equation (7) can be written in 
terms of quantities that can be calculated as 


u+t 


U+ Vm 


r++ 


y 


m 


(7a) 


where 


U + = F(y + ) = t £>)-»(£ ft) « 

The function F is obtained from the relation between u + and y + in 
figure 1. The parameter r ++ is a type of Reynolds number and is as- 
signed an arbitrary value. In order to calculate y + , a force balance 
is written on the element AA in the preceding sketch. The forces on 
the element are the shear force acting on Al and the pressure forces 
acting on the faces of the element. There are no shear forces acting 
on the velocity-gradient lines because the normal velocity derivatives 
are zero along those lines. Writing a force balance on the element 
gives 


A. 


= _AA(<3£) 

T 0 A l \dx / 


( 9 ) 
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where the pressure gradient dp/dx is uniform over the annulus because 
the flow is fully developed. Substituting equation (9) into the defini- 
tion of y+ results in 

4 - r++ t do) 

where AA, Al , y , and r are measured from a figure similar to the 

preceding sketch and r"*"*" is a parameter. If, for a specified value of 
r++, we consider a given line normal to the surface, the value of y/y m 

for a given u*'*' can be calculated by using equations (10) , (7a) , and 
(8) in that order. 

By carrying out the calculation for various lines normal to the 
surface, lines of constant velocity (constant u ++ ) can be obtained. 

As mentioned previously, new and more accurate velocity-gradient lines 
are next drawn so as to intersect the constant-velocity lines at right 
angles. The calculation is then repeated with the new velocity -gradient 
lines . 


Calculated Velocity Distributions 

Several velocity distributions for flow in equilateral triangular 
and in square passages, calculated by the method described, are shown 
in figure 3. The calculations were carried out for severed. Reynolds 
numbers, but the effect of Reynolds number was found to be negligible. 
In all cases the constant-velocity lines (solid lines) are essentially 
normal to the velocity-gradient lines (dashed lines) as required. Com- 
parison of these constant -velocity curves and the experimental curves 
of Nikuradse (ref. l) indicates that both have the same general shape, 
although the experimental curves seem to approach the corners more 
closely than the analytical ones. This might indicate the presence of 
secondary flows in the corners, as originally suggested by Frandtl, 
but it appears that no definite conclusions can be drawn on the basis 
of the available experimental data. 


FRICTION FACTOR AND REYNOLDS NUMBER 


With the velocity distributions for the passage known, friction 
factors and Reynolds numbers can be calculated by integrating the distri- 
butions to obtain bulk or average velocities. The bulk velocity between 
two adjacent straight lines normal to the wall (dashed lines in preceding 
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sketch) is first obtained. This bulk velocity varies with position on 
the tube and is given by 


"b = 



u dA' 


AA' 


( 11 ) 


where AA’ is the area between two adjacent straight lines normal to the 
wall. Equation (ll) can be written in dimensionless form as 



The average velocity for the whole passage is 


( 12 ) 



dA 


(13) 


where A is the total area of the passage." The variation of 

u++/u++ = u. /u, with dimensionless distance from the corner for a 

h ' p,av D' D,a,v 

low and a high Reynolds number is shown in figure 4. 


For later comparison with heat transfer, the ratio of the local shear 
stress to the average shear stress can be obtained from equation (9), and a 
corresponding equation for the whole passage, as 

T 0 _ AA _r_ 

T 0,av A AZ 

Curves for Tq/tq av against dimensionless distance from the corner 

are presented in figure 5. The curves indicate that the shear stress is 
lowest in the region near the corner and goes to zero at the corner. The 
slopes of the Tq/tq av curves near the corners are not as steep as those 

of Nikuradsej this is probably due to secondary flow in the corners, which 
is not considered in this analysis. The value of Reynolds number has a 
small effect on the shape of the curves. The values for the Reynolds 
numbers were calculated from 


(14) 


Re 


-H* r++ 

av 






(15) 


which follows directly from the definition of Reynolds number. The hy- 
draulic diameter D e is defined in the usual way as four times the flow 
area over the wetted perimeter. 

The fact that the shear stress should decrease in the region near the 
corner can be seen directly from equation (9), which indicates that the 
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shear stress is proportional to AA/AZ. But AA ~ AZy m , and, hence, the 

shear stress is approximately proportional to y m , which decreases as the 

2 

corner is approached. 

The friction factor for the passage, Based on the pressure gradient, 
is defined "by 


or, in terms of dimensionless groups, 

* ^-r 

2(r/ D e )UbV 

Figures 6(a) and (b) show friction factors based on the pressure 
gradient plotted against Reynolds number for the equilateral triangular 
and the square passages, respectively. Data from reference 10 are in- 
cluded in the figure. The data for the square passage are in good agr ee - 
ment with the predicted curve; those for the triangular passage are 
about 10 percent above the predicted curve but are still below the 
circular-tube line taken from reference 7. 

Figures 6(c) and (d) show analytical and experimental friction 
factors for axial flow between rods . These experimental values (from 
ref . ll) were not available at the time the analysis for flow between 
tubes was presented (ref. 3) and so are included herein. The original 
results in reference 11 have since been recomputed by the authors of that 
report, and the corrected data are shown in figures 6(c) and (d) . Th^ 
agreement between analysis and experiment is reasonably good. It is 
possible that the higher values of the friction factors at the higher 
Reynolds numbers in the case of the square array are caused by secondary 
flow. It was found in reference 12 that secondary flows occurred only 
at the higher Reynolds numbers. 




D e dp/dx 

2 p uj? 

H D ; av 


( 16 ) 

(17) 


WALL HEAT-TRANSFER 013151181111011 

With the velocity distributions previously obtained (and with certain 
assumptions) , the fully developed heat-transfer distributions are calcu- 
lated herein. From these values, the wall temperature distributions 
will be calculated in the next section. 

The heat added between two straight lines normal to the wall per 
unit length of tube is q^Z (see sketch on p. 5). It is assumed herein 

‘‘It is perhaps surprising that the curves in figure 5 for the triangu- 
lar passage, as well as those in figure 3(a), have nearly the same shape as 
corresponding curves calculated for l amin ar flow. 
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that all this heat is used in heating the fluid element between the 
straight lines. This is a good assumption when the peripheral, tempera- 
ture variation is small compared with the variation of temperature 
across the passage. 


Making a heat balance on the element of fluid between two straight 
lines normal to the surface and using the foregoing assumption result in' 3 


q^l = p^Cp 



AA* 


For the whole passage cross section. 


( ^0,av r - P u b,av c p 



(18) 


(19) 


It is shown in appendix B that, for the fully developed case, 
dt^/dx = dt^ av /dx when the heat transfer per unit area at* a given cir- 
cumferential location does not vary with _x. Division of equation (18) 
by equation (19) and conversion to dimensionless form give 


go 

^0,av 


n TT 

u b AA r 
, .++ A AZ 
D,av 


( 20 ) 


Equation (20) also gives the ratio of local Husselt number to av- 
erage Nusselt number for the case of uniform circumferential wadi tem- 
perature, when the heat-transfer coefficients sire defined in the following 
way: 


_g0 % 

t 0 - *b ^0,av " ^av 


( 21 ) 


and 


, _ 3o,av 

aV = ^0,av - *b,av 


( 22 ) 


or 


d- . ja. . Al_ (23) 

h av Nu av lo,av 

■z 

As an alternative assumption, an element bounded by velocity-gradient 
lines rather than by straight lines was used. No net heat transfer across 
the velocity-gradient lines was assumed. This assumption was found to 
give essentially the same distribution of as did the assumption. _used 

in the text. 
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Equation (23) is true only for the case of uniform circumferential wall 
temperature, "because for any other temperature distribution the tempera- 
ture differences in equations (2l) and (22) will not cancel. 

Values of qo/lO,av (^Aav for ^ fozm circumferential wall temper- 
ature) are shown as a function of dimensionless distance frcm the corner 
in figure 7. As was the case with the shear stress (fig. 5), the heat 
transfer decreases as the corner is approached and goes to zero at the 
comer. However, comparison of figures 5 and 7 indicates that the heat 
transfer approaches zero much more rapidly than does the shear stress. 
Comparison of these figures indicates, therefore, that the assumption 
of similarity of the shear stress and heat-transfer distributions, which 
is made in reference 1, might be expected to give heat-transfer coeffi- 
cients in the vicinity of a comer that are too high. Actually, that 
assumption gives results that are in reasonable agreement with experiment 
(ref. 10), possibly because of the compensating effects of secondary 
flows and peripheral heat transfer in the fluid. 

The reason that the heat transfer decreases more rapidly than the 
shear stress as the corner is approached can easily be seen by comparing 
equations (9) and (18). Both equations contain AA/AZ, or AA'/AZ, which 
causes the shear stress or heat transfer to decrease near the corner. 
However, the heat-transfer equation (eq. (18)) contains, in addition, 
the local bulk velocity u^, which also decreases j therefore, the heat 

transfer decreases more rapidly than the shear stress as the corner is 
approached. 


WALL TEMPERATURE DISTRIBUTION 

It is assumed in this section that the passage walls are thin and 
that heat is transferred uniformly to the outside surface of the passages. 
Uniform heat sources may also be present in the wall. Because of tangential 
conduction around the wall, the heat transfer through the inner surface 
will not be uniform and will have the distributions obtained in the pre- 
ceding section on wall heat transfer. The heat transfer per unit area 
through the outer surface plus the heat generated per unit volume times 
the wall thickness is equal to q^ , the average heat transfer per 

unit area through the inner surface to the fluid. 

In order to obtain the temperature distribution around the passage, 
a heat balance is first made on an element of wall of circumferential 
length dZ. This heat balance gives 


v 

^0,av *() dZ 


(24) 
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where b is the thickness of the wall and q^ is the tangential heat 

conduction per unit area. Equation (24) can he written in integral form 
as 



where q^. is zero at the corner because the wall temperature distribution 
is symmetric about that point (temperature gradient is zero) . However, 

dt 0 

= “ k t dZ ( 26 ) 

Substituting equation (26) in (25) and integrating again result in 



where t n ^ is the wall temperature at the corner. By using the values 
of Qq/<3q aT obtained in the preceding section, the difference between 

the maximum wall temperature and the wall temperature at any point for 
a given heat flux can be calculated from equation (27). A dimensionless 
parameter containing the difference between the maximum and the average 
wall temperatures can be obtained by integrating equation (27); that is, 



The relation for k^b(tQ - tQ^ av )/qQ^ av r 2 can be obtained by subtracting 
equation (27) from (28). 

The results for the wall temperature__distributions are presented in 
figure 8, where k t b(tQ - tQ^ ay )/qQ^ av r 2 is plotted against dimension- 
less distance from the wall. As mentioned previously, the quantity qQ av 
in the wall temperature parameter is equal to the sum of the uniform heat 
flux at the outer wall of the tube and the heat generated per unit volume 
times the tube wall thickness. Figure 8 indicates that t 0 ^ raax - t 0 ^ av , 

which is the quantity of greatest practical interest, is directly propor- 
tional to qg^ av . 
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Experimental wall temperature differences from reference 10 were 12 
to 60 percent lower than those in figure 8. This difference might "be 
caused "by secondary flows or "by peripheral conduction and turbulent 
transport in the fluid, both of which were neglected in the analysis, or 
by nonuniform heat generation in the test- section walls. 


AVERAGE HEAT-TRANSFER COEFFICIENTS 

All the results obtained thus far (wall heat-transfer distributions, 
wall temperature distributions, etc.) were independent of Brandt 1 number, 
inasmuch as it was not necessary to use the generalized temperature dis- 
tribution in figure 2. However, it is necessary to use that distribution 
for obtaining the difference between wall and bulls temperatures that 
corresponds to a given heat flow,' that is, for obtaining the heat-transfer 
coefficient . 

The average heat-transfer coefficient for the inner wall of the 
passage is defined by 


where 


and 



(29) 


(30) 


(31) 


(32) 


r T T ht= average Nusselt number corresponding to the average heat-transfer 
coefficient (eq. (29)) can be calculated from the following equation, which 
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can be verified by substituting the definitions of the various quantities: 
1 r / D e I % 


where 


and 


1 5 = •»"&)*&) 


^ U av r+^Tr ^o,&v 

r kr I k t b(t 0 - tg^av) / ^ 

Dgkt'b I « -2 


'0 


^av rt 


(^-) a (±) 

\ u b , av/ \^0/ 


(33) 


V - 


r 1 

j.++ u++ / A \ 

^ 4 W 


(34) 


t++ 


t + t ++ y m /r 


(35) 


m 


As in the case of velocity distributions, the values of y + are measured 
along lines that are normal to the wall and extend to the line of maxi- 
mum velocities. The values of t + are obtained from figure 2. 

Average Nusselt numbers calculated by equation (33) axe plotted 
against Reynolds number for various values of the parameter kr/kj/b and 

of Prandtl number in figure 9. The parameter kr/ktb is zero for 

uniform wall temperature and increases as the variation of wall tempera- 
ture increases. The wall temperature distributions were obtained from 
figure 8. The Nusselt numbers decrease as kr/k^b increases or as 

tube conductivity or thickness decreases^ that is, the Nusselt numbers 
decrease as the wall temperature variation increases. These trends are 
qualitatively similar to those obtained for laminar flow in reference 13, 
where it was found that the Nusselt number for a rectangular duct with 
uniform peripheral heat flux was lower than that for uniform wall 
temperature . 

Experimental data from reference 10 are included in figures 9(d) 
and (h) . In general, the data are in reasonably good agreement with 
the predicted values. It Is possible that the slight departure from 
the predicted curves at the higher Reynolds numbers is caused by secondary 
flow (ref. 12). 
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SUMMARY OF RESULTS 

The following results were obtained, from the analytical investigation 
of fully developed axial turbulent flow in noncircular passages with heat 
transfer : 

1. The velocities and shear stresses in the region near the corner 
were lower than the average values and went to zero at the corner. 

2. When the passages were heated, the heat transfer to the fluid in 
the region near the corner was lower than the average value and went to 
zero at the corner. 

3. The friction factors for the noncircular passages were somewhat 
lower than those for a circular tube. 

4. When uniform heat sources in the passage wall and uniform heat 
transfer at the surface were assumed to occur, the difference between 
the maximum and average wall temperatures was directly proportional to 
the heat flux. 

5. The average Nusselt numbers for the noncircular passages were 
somewhat lower than those for a circular tube. The average Nusselt 
number was also found to be a function of wall temperature distribution 
and of Prandtl number. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 21, 1958 
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APPENDIX A 


SYMBOIS AND DIMENSIONLESS QUANTITIES 


total area of section of passage 
wall thickness of passage 

specific heat of fluid at constant pressure 


equivalent diameter 


friction factor. 5 

2 e“b,av 

<30 

local heat-transfer coefficient, x r — 

t 0 " %,i 


average heat -transfer coefficient. 


thermal conductivity of fluid 


' fc 0 ,av " "^h,! 


thermal, conductivity of material of passage wall 


wall temperature distribution 


peripheral distance along the wall from corner 
Nusselt number based on local heat transfer, hD e /k 
Nusselt number based on average heat transfer, ■^av^e/^ 

constant , 0 . 124 
Prandtl number, Cp|J./k 
static pressure 
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rate of heat transfer per unit area 

rate of heat transfer from wall of passage per unit area 

average rate of heat transfer from wall of passage per 
unit area 

Reynolds number, pu^Dg/p 

distance from corner to midpoint of wall 


l-rfdp/dx) 


parameter, 


temperature of fluid at a point 

(*0 " t ) c p T 0 


temperature parameter, 


'Jq y/r/p 


temperature parameter, 


(to - t)cpp 


l-rfdp/dx) 


arhitrary temperature at a given cross section 


local "bulk temperature of fluid at given distance from 
corner at cross section of flow passage 


hulk-temperature parameter 


•C ’»"$*(&) 


wall temperature 

maximum wall temperature of passage 


time -average velocity parallel to wall at a point 


velocity parameter, u/ V^Tq/p 


velocity parameter. 


-r (dp/ dx) 
P 
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% 



X 

y 

y + 

ym 

y+ 

m 

a 

e 

®h 

x 

n 

v 

P 

T 

Subscripts : 

av 

t 


local bulk velocity at given distance from corner at 
cross section of flow passage 


bulk -velocity parameter, 

f A 

bulk-velocity parameter, i I dA 



axial distance along passage 
normal distance from wall 


wall distance parameter, 


V T o/p 

p/p 


y 


value of y at u = u m 

y+ evaluated at line of maximum velocities 


ratio of eddy diffusivity for beat transfer to eddy 
diffusivity for momentum transfer, %/ e 

coefficient of eddy diffusivity for momentum 
coefficient of eddy diffusivity for heat 

Karman constant, 0.36 
absolute viscosity of fluid 
kinematic viscosity 
mass density of fluid 
shear stress in fluid 


0 


average 

tangential 

pertaining to a wall 
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APPENDIX B 

CONDITION FOE FOLLY DEVELOPED BEAT TRANSFER 
For fully developed heat transfer and constant fluid properties , 

t 0 - t a n a 

where 1%, an arbitrary heat-transfer coefficient, is independent of x. 
If h a were not independent of x at a great distance from the entrance 
(cyclic variations of h a excluded) , the absolute value of h a would 
become arbitrarily large as x increased, so that for finite temperature 
differences the absolute value of < 1 q would become arbitrarily large. 

The following equations are special cases of equation (Bl) : 



Qo 

t 0 " *b,av 


(B2) 


90 
^O " 


= h 


(B3) 


For the case where the wall heat transfer per unit area is independent 
of x (but not of circumferential position), equations (B2) and (B3) 
can be differentiated to give the following results: 

dt 0 = dt b,av 
dx dx 

«0_«b 
dx dx 

or 

dt b _ dt b,av 
dx dx 

The quantity dt^/dx is therefore independent of circumferential 
position when the heat transfer Is independent of x. 
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Figure 1. - Generalized velocity distribution for adiabatic turbulent flow. (Vertical dashed line is dividing line 
■between eqp« ) 
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(a) Triangular passage. 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Dimensionless distance from corner, i/r 

(b) Square passage. 


Figure 4. - Predicted variation of bulk velocity parameter with distance from corner. 
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Figure 5. - Predicted variation of shear stress on wall with distance from corner. Reynolds number. 
24,000 or 900,000. 
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Friction factor 



(a) Triangular passage. 



Reynolds number, Re 
(b) Square passage. 


Figure 6. - Predicted friction factors based on static-pressure gradient as a 
function of Reynolds number. 
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Figure 6. - Concluded. Predicted friction factors based on static-pressure 
gradient as a function of Reynolds number. 




Heat- transfer parameter, qg/q 0 


30 


WACA TN 4384 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Dimensionless distance from corner, i/r 


Figure 7. - Predicted variation of heat transfer from passage wall with distance from corner. 
Reynolds number, 24,000 or 900,000. * T " 


4949 




Vail temperature distribution 
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Figure 8. - Predicted variation of wall temperature of passage with distance from corner. Reynolds 

number, 24,000 or 900,000. iH 
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(c) 'triangular passage, ^r- 
distribution) . *t D 


0.010 (nonuniform peripheral wall tecgp. 


Figure 9- - Continued. Predicted airerage Nub salt number as a funct 


(d) Triangular passage, * 0.025 (non 


uniform peripheral wall temp, 
distribution) . 


of Reynolds number and Prandtl number. 
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(h) Square passage, - 0,015 (nonunlforra peripheral 
wall temp, distribution). 


a funotlon of Reynolds number and Prandtl number. 
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